Parametric generation of quadrature squeezing of mirrors in cavity optomechanics 
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We propose a method to generate quadrature squeezed states of a moving mirror in a Fabry- 
Perot cavity. This is achieved by exploiting the fact that when the cavity is driven by an external 
field with a large detuning, the moving mirror behaves as a parametric oscillator. We show that 
parametric resonance can be reached approximately by modulating the driving field amplitude at a 
frequency matching the frequency shift of the mirror. The parametric resonance leads to an efficient 
generation of squeezing, which is limited by the thermal noise of the environment. 
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I. INTRODUCTION 

Cavity optomechanics [lJ-Q , as an interaction inter- 
face between a cavity field and a moving mirror, is an 
exciting research area for exploring quantum behavior 
in macroscopic systems as well as applications in quan- 
tum information processing. With the recent advances 
of cooling techniques in optomechanical systems |5l— ll 2jj , 
it is becoming possible to overcome thermal noise and 
study quantum state engineering of mechanical mirrors. 
Indeed, recent studies have already shown that various 
kinds of non-classical states can be generated by optome- 
chanical coupling. These include quan tum superposition 
states [HI, H], entangled states I15l-fl9ll. and squeezed 



states of light [2fJ-|23| and mirrors 0128 1. 

Specifically, achieving squeezed states in mechanical 
oscillators (mirrors) is an important goal because of the 
applications in ultrahigh precision measurements such 
as the detection of gravitational waves mum. Several 
schemes have been proposed to create quantum squeezing 
of the moving mirror in cavity optomechanics. For exam- 
ple, squeezing can be transferred from a squeezed light 
driving the cavity to the mirror ;25j, and recently Mari 
and Eisert have shown that squeezing can be generated 
directly by a periodically modulated driving field [26]. 

We note that a basic mechanism for creating quadra- 
ture squeezing is to introduce a parametric coupling for 
the motional degree of freedom of the mirror. In particu- 
lar, efficient squeezing can be achieved at the parametric 
resonance, such that the Hamiltonian in the interaction 
picture takes the form Hj cx b 2 + b^ 2 [where b and w are 
operators of the oscillator in Eq. (JTJ] and the correspond- 
ing evolution operator is a squeezed operator. Therefore 
an interesting question is how the parametric resonance 
can be reached in cavity optomechanical systems. One 
of the difficulties here is the dynamical shift of the me- 
chanical resonance frequency due to the optomehanical 
coupling, which is sensitive to the intensity of the cavity 
field. In this paper we show that in the large detuning 
limit, the frequency shift can be compensated by modu- 
lating field amplitude at a suitable frequency, and hence 
parametric resonance can be reached approximately. We 
will present an explicit form of the driving amplitude, 



and analyze the time development of squeezing in the 
presence of thermal noise. 



II. MODEL 

The system under consideration is an optical cavity 
formed by a fixed mirror and a moving mirror connected 
with a spring (Fig.[T]). We consider a single-mode field in 
the cavity and model the moving mirror as a harmonic 
oscillator. The Hamiltonian of the system reads 



Hs = huj c a} a + Hu) m b'b — hga) a{b^ + b) 
+htt{t)e- iu)dt a i + hfl*{t)e lUdt a, 



(1) 



where (b*) and a (b) are the creation and annihilation 
operators associated with the single-mode cavity field 
(mirror) with frequency u c (w m ). Assuming m e s is the 
effective mass of the mirror, then the position and mo- 
mentum operators of the mirror are x — x zp { (p* + b) and 
p = irn ffW m a; Z pf(6 t - b), where x zp{ = ^h/(2m cS u> m ) is 
the zero-point fluctuation of the mirror's position. The 
third term in Eq. (fT]) describes a radiation pressure cou- 
pling with the coupling strength g — uj c x zp i/L, where L 
is the rest length of the cavity. In addition, the cavity is 
driven by an external field with a main frequency uid and 
the time- varying amplitude Q(t). 

To include damping in our model, we follow the stan- 
dard approach by coupling the system with oscillator 
baths such that the quantum Langevin equations (in a 
rotating frame with frequency w<j) for the operators a and 
b are given by 

a = —iA c a + iga(w + b) — iCl(t) — + a in , (2a) 



b = -iu m b + iga'a- ^-b + b in , 



(2b) 



with the detuning A c = lj c — Ud and the cavity (mirror) 
decay rate j c { r y m )- Under the assumption of Marko- 
vian baths, the noise operators a.;„ and bi n have zero 
mean values and they are characterized by the correlation 
functions (a in (t)a\ n {t')) = J c 5(t-t'), (a\ n {t)a ln {t')) = 0, 
(bin{t)bUt')) = lm {n m + l)8{t-t'), and (bUt)b in {t')) = 
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FIG. 1: (Color online) Schematic diagram of the cavity op- 
tomechanical system. An externally driven Fabry-Perot cav- 
ity is formed by a fixed end mirror and a harmonically bound 
end mirror. 



l m n m 5(t - 1'), where n m = {exp[hu m / (k B T m )} - 1} 1 is 
thermal excitation number of the mirror's bath at tem- 
perature T m and fcs is the Boltzmann constant. Here we 
have assumed that the bath coupled to the cavity field 
is effectively a vacuum, and the rotating-wave approxi- 
mation has been employed to describe the system-bath 
interaction 133. 133I1. 



III. LINEARIZED SYSTEM: FORMAL 
SOLUTION OF FLUCTUATIONS 

Next we write a = {a} + Sa and b = (b) + Sb such 
that the fluctuations about the expectation values are 
described by operators Sa and Sb. Assuming the fluctua- 
tions are sufficiently small, then we may linearize Eq. ([2]) 
to obtain the equation of motion for Sa and Sb: 



Sa = -iA(t)Sa + ig{a{t)){8tf + Sb) 



^-Sa + a in , (3a) 

Sb = -ioj m Sb + ig[(a\t))Sa + (a{t))5a?} - ^5b + b in , 

(3b) 

where A(i) = A c - g[{b(t)) + (& f (t))]. The expec- 
tation values (a(t)) and (b(t)) are governed by equa- 
tions of motion: (a) = — [iA(t) + -j](a) — iil(t) and 

(6) =-(^ m +~ J f)(b)+ig\(a)\ 2 . 

For convenience, we introduce the quadrature opera- 
tors by 8X s = aib = (Ss^ + <5s)/v2 and SY s = a , b = i(8s^ — 
6s)/y/2. Then Eq. (3) can be concisely expressed as 



v(t) = M(t)v(t) + N(t) 
where v = (SX a ,SY a , 5X b , SY b ) T , and M is 



(4) 



M(f) 



A(t) -V2g(Y a (t)) 
V2g(X a (t)) 



V2g(X a (t)) V2g(Y a (t)) 



_2£ 

-A(t) 




( 5 ) 

with (X s = a>b (t)) = [< S t(f)) + ( s ( t ))]/V2and (r s=0 , 6 (t)> = 
*[(s t (*)) - (s(*))]/V^- The noise vector in Eq. @ is 



defined by N = (X l a n ,Y™,X b in ,Y b in ) T , with X*IL 0i6 = 
(4 + s m )/V2 and y;« ai6 = i(4 - Sm )/V2. 

Equation Q is a first-order linear inhomogeneous dif- 
ferential equation with variable coefficients. Its formal 
solution is 

v(i) = G(t)v(0) + G(t) / G- 1 (r)N(r)dr, (6) 



where the matrix G(£) satisfies G(i) = M(t)G(i) and 
the initial condition G(0) = / (/ is the identity ma- 
trix). In the present system, interesting quantities are 
the quadrature fluctuations of the cavity and the mir- 
ror. Hence, we define a covariance matrix R(i) by the 
elements Rjj'(t) = (vj (t)vj/ (*)) for = 1,2,3,4. Obvi- 
ously, the four diagonal elements of R(i) are the expecta- 
tion values of the square of the four quadrature operators 
of the system. They are Rn(t) = (8X%(t)), R22O) = 
(5Y a 2 (t)), R 33 (t) = (SX 2 (t)), and R 44 (i) = (SY b 2 (t)). 
For the mirror's rotating quadrature operator X b (9,t) = 
cos6 , ATh(t) + sin0Y(,(i), the corresponding variance is 
given by (SX 2 {6,t)) = cos 2 6»R 33 (t) + sin 2 6»R 4 4(t) + 
i S in20[R 34 («)+R 43 (f)]. Since [X b (6,t),X b (6+ir/2,t)] = 
i, quadrature squeezing occurs when (SX 2 (9 7 t)) < 1/2. 

To test the dynamical quadrature squeezing, we need 
to determine the covariance matrix R(t), which has the 
formal expression: 

R(i) = G(i)R(0)G T (i) + G(t)Z(t)G T (t). (7) 

where Z(t) is defined by 

Z(t)= f f G-\T)C{T,T')[G-\r')] T dTdT'. (8) 
Jo Jo 

Here C(r, r') is the two-time noise operator correla- 
tion matrix defined by the elements: C„ n '(r, r') = 
(N„(r)N n /(r')) for n, n' = 1,2,3,4. For Markovian 
baths, we have C(r, r') = C5(r — r'), where the con- 
stant matrix C is given by 
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(9) 



GENERATION OF QUADRATURE 
SQUEEZING 



Having obtained the formal equations for the evolu- 
tion of quadrature fluctuations of the mirror, we now ask 
how the external driving amplitude Cl(t) can be chosen 
to generate a large degree of quadrature squeezing of the 
mirror. We approach the problem by considering the 
large detuning regime (A c S> u> m ) so that by adiabatic 
approximation we have 



Sa 



A c - i lc /2 



[a(t))(6tf + 8b) + ft, 



(10) 
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FIG. 2: (Color online) Time evolution of (5X 6 2 (tt/4, t)) at var- 
ious temperatures T m . From the bottom up, the three curves 
correspond to kBT m /(frw m ) = 0, 20, and 50, respectively. The 
standard quantum limit 1/2 is indicated by dashed black line. 
The parameters are given in the text. 

with F in = r o t a m (i')e (iAc+7 <= /2)(t '" t) dt'. Here, we have 
also assumed A c ^> g(Xf > {t)) and hence A(t) sa A c . Cor- 
respondingly, the equation of motion (|3bj) for Sb becomes 

Sb = -iu m 5b + i V \(a(t))\ 2 (Stf + Sb) - ^Sb + F[ n , (11) 

where rj — ^2^.^2/4 an d the noise operator consists of 
two parts F' m = F t a n + b m . The part f\° = ig (a* (t ) ) F m + 
ig{a(t))F- n comes indirectly from the cavity's bath and 
depends on the mean field solution, while the second part 
bin comes directly from the mirror's bath. 

Next we observe that if the external driving amplitude 
is chosen as 

fi(t) = fl sin -&)*], (12) 

with Q being a constant and £ = g 2 ^ 2 l A c /(A 2 .+7 2 /4) 2 , 
then by the adiabatic solution (a(t)) ps — f2(t)/(A c — 
ijc/2) and the assumption tu m 3> £bj Eq. (|lip can be 
approximated by 

SB = -i^-SB^ - ^SB + F' ln e l ^-^\ (13) 

where SB = Sbe l ^ Um ~^ t is defined. In deriving Eq. ([15]). 
we have made use of a rotating wave approximation such 
that counter-rotating terms with the rapidly oscillating 
phase factors e ±2i(a>m-fo)t and e ±4i(w m -£ )t nave been 

dropped. 

We notice that Eq. (jT3l precisely corresponds to the 
equation of motion of a damped parametric oscillator at 
resonance. If damping can be ignored, a mirror initially 
prepared in the ground state would display exponen- 
tial squeezing as time increases: (SX 2 (tt/A, t)) — ^e - ^ *. 
Such an efficient squeezing can be understood by inspect- 
ing Eq. (fTTj) in which our choice of f2(t) matches the aver- 
age value of the shifted resonance frequency of the mirror 



w m — 7?|( a (£))| , and therefore the parametric resonance 
can be reached approximately. Note that £0 is the aver- 
age value of such a frequency shift and it also plays the 
role of an effective strength of the parametric process. 

However, it should be noted that for practical pur- 
poses, £0 in Eq. ([T3| . which decreases as A~ 3 , has to 
be strong enough to overcome noises of the baths, i.e., 
the detuning A c cannot be arbitrarily large. For real- 
istic choices of A c , the quality of squeezing has to be 
examined in the presence of noise without making use 
of the adiabatic approximation. To this end, we em- 
ploy the linear formalism above and solve directly the 
covariance matrix in Eq. (J7J numerically. For simplic- 
ity, we assume that the system is initially prepared in 
its ground state |0) c ® |0) m through a state prepara- 
tion process. Such an initial state may be achievable in 
future experiments based on ground-state cooling tech- 
niques. In addition, we consider the following systems 
parameters: u m — 2tt x 1 MHz, A c = 2tt x 10 MHz, 
7m = 2tt x 100 Hz, 7c = 2?r x 100 kHz, fl « 31.6 GHz 
and g = 2ir x 100 Hz, which are realistic under current ex- 
perimental conditions [34|,[35j]. In Fig. [3] we plot the time- 
dependence of quadrature variance of the mirror at vari- 
ous temperatures based on the form of O(i) in Eq. (|12p . 
the evidence of squeezing is clearly shown at sufficiently 
low temperatures. In fact, for not too large A c = 10aj m 
chosen in Fig. [3J our exact numerical results agree well 
with the adiabatic approximation. 

If the temperature of the mirror's bath is higher than a 
critical value, X^, then there will no longer be squeezing 
in the mirror (Fig. (2[ blue line) . A rough estimation of 
the damping effect can be made by considering that the 
noise is mainly from the mirror's bath, so that 

(5X?(ir/i,t)) a l e -(7 m +fr)* + 7m("m + l/2) / e _ (7m+fo)t \ _ 
2 (7 ?n + £0) v ' 

(14) 

Therefore squeezing occurs if the thermal excitation num- 
ber n m is below a critical number, n c m = For the 
parameters used in Fig. [21 our estimation gives 
4.8 mK or kBT„\/tko m ~ 50.5, and this agrees with the 
numerical value 50 shown in Fig. [5J We remark that in 
deriving Eq. (14), the effect of noise aj n has been ne- 
glected. This can be justified by a lengthly calculation 
which shows that (SXf) due to the cavity field noise is of 
the order of £o(7m + £0 + 7c)/[4A c (7,„ + £ )] in the long 
time limit, and hence it can be made small compared with 
the contribution from the thermal bath of the mirror by 
a large detuning. 



V. CONCLUSION 

To conclude, we have presented a method to generate 
quadrature squeezing of a mirror in cavity optomechan- 
ics. Specifically, we have shown that in the large detun- 
ing regime with A c > uj m > £ and A c > g\(X b (t))\, 
the driving field of the form f2(t) given in Eq. (12) can 



4 



generate squeezing dynamically 36]. The squeezing is 
supported by direct numerical calculations for realistic 
parameters. We should point out that our scheme is dif- 
ferent from that in Ref. [26} because the large detuning 
regime considered here enable us to eliminate the cav- 
ity field and formally map the mirror to a parametric 
oscillator. In addition, parametric resonance can be fine 
tuned by our driving field Cl(t) so that the frequency shift 
of the mirror due to coupling to the cavity field can be 
compensated approximately. 
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